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Solid "^He and the Supersolid Phase: 
from Theoretical Speculation to the Discovery of a New State of Matter? 
A Review of the Past and Present Status of Research 

Davide Emilio Galli* and Luciano Reatto'^ 

Dipartimento di Fisica, Universitd degli Studi di Milano, Via Celoria 16, 20133 Milano, Italy 

The possibility of a supersolid state of matter, i.e., a crystalline solid exhibiting superfluid 
properties, first appeared in theoretical studies about forty years ago. After a long period of 
little interest due to the lack of experimental evidence, it has attracted strong experimental 
and theoretical attention in the last few years since Kim and Chan (Penn State, USA) reported 
evidence for nonclassical rotational inertia effects, a typical signature of superfluidity, in samples 
of solid *He. Since this "first observation", other experimental groups have observed such 
effects in the response to the rotation of samples of crystalline helium, and it has become clear 
that the response of the solid is extremely sensitive to growth conditions, annealing processes, 
and ^He impurities. A peak in the specific heat in the same range of temperatures has been 
reported as well as anomalies in the elastic behaviour of solid ''He with a strong resemblance to 
the phenomena revealed by torsional oscillator experiments. Very recently, the observation of 
unusual mass transport in hep solid ''He has also been reported, suggesting superflow. From the 
theoretical point of view, powerful simulation methods have been used to study solid ^He, but 
the interpretation of the data is still rather difficult; dealing with the question of supersolidity 
means that one has to face not only the problem of the coexistence of quantum coherence 
phenomena and crystalline order, exploring the realm of spontaneous symmetry breaking and 
quantum field theory, but also the problem of the role of disorder, i.e., how defects, such as 
vacancies, impurities, dislocations, and grain boundaries, participate in the phase transition 
mechanism. 

KEYWORDS: supersolid, quantum solid, solid ''He, nonclassical rotational inertia, superfluidity, Bose- 
Einstein condensation, off-diagonal long-range order, zero-point motion, disorder, vacancies 
zero-point vacancies, interstitials, vacancy-interstitial pairs, dislocations, grain boundaries 
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1. Introduction 

Should we see atoms flowing without resistance in and 
out of a solid made of the same atoms, this would cer- 
tainly be counterintuitive because the same entities have 
to simultaneously give the rigidity intrinsic to a solid and 
provide the moving particles. Note the difference between 
this system and the case of a superconductor, where the 
flowing charges are the electrons and the solidity of the 
system is given by the ions. Still, this possibility is what 
should happen in a supersolid, a state of matter conjec- 
tured to exist almost forty years ago,''^-' and solid *He 
appeared to be the most likely candidate. A plethora 
of counterintuitive behaviours have their origin in quan- 
tum mechanics and in the indistinguishability of identical 
particles, and supersolidity also has its roots in quantum 
mechanics. In the case of identical particles with integer 
spin, such as *He, which has zero total spin, the exten- 
sive occupation of a single quantum state, i.e., the phe- 
nomenon of Bose-Einstein condensation (BEG) , is at the 
heart of (in the sense that BEG is proven to imply'^^) 
superfluidity. BEG has been recognized in the superfluid 
phase of liquid ^He, and more recently it has been di- 
rectly observed in the metastable phases of low-density 
ultracold trapped atomic systems. Typical characteris- 
tics of the superflow in such phases^-' are nonclassical 
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rotational inertia (NGRI) and also the absence of dissi- 
pation for flow velocities under the critical velocity as 
well as the presence of vortices with a velocity flow field 
with quantized circulation in units of h/m, where h is 
Planck's constant and m is the mass of the particles. 

As already mentioned, a solid phase with such super- 
flow properties seems paradoxical at flrst sight. However, 
from a theoretical point of view, nothing prohibits such 
an occurrence. In fact, model quantum systems exist in 
which the existence of a supersolid phase has been proved 
(an example inherent to the topic of this review is pre- 
sented in ref. 5) or found (for quantum lattice models 
see, for example, refs. 6-8). A supersolid phase would 
correspond to a phase in which two kinds of order are 
simultaneously present: crystalline long-range order and 
off-diagonal long-range order^) (ODLRO). Thus, in such 
a state, both order in real space (crystalline order) and 
order in momentum space (corresponding to the off- 
diagonal order) should be simultaneously present. These 
two different kinds of order originate from two differ- 
ent broken symmetries: the translational symmetry, as 
observed in any crystalline solid (directly detectable in 
the Bragg scattering), and the gauge symmetry, whose 
breakage allows for the phase coherence throughout the 
system. Our discussion is devoted to solid '*He, in which 
the translational symmetry is spontaneously broken and 
only the interatomic interaction plays a role in determin- 
ing the reference lattice. We are not going to discuss the 
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cases of adsorbed phases of helium, of bulk systems of 
different species such as hydrogen, of quantum gases in 
optical lattices, or of lattice models. Note that lattice 
models are qualitatively different from the case of ^He, 
because in such models the Hamiltonian is invariant only 
for lattice translations so that the spatial broken symme- 
try refers to a discrete group, not to a continuous group 
such as in a quantum solid. 

In this review we have not tried to produce an exhaus- 
tive report of the research contributions to this (presently 
hot) topic, which are rapidly growing in number. We 
suggest that the reader also refers to other reviews^°~^^^ 
to obtain a more comprehensive view of what has been 
achieved and what remains to be achieved in the theo- 
retical and experimental investigation of a possible su- 
persolid ^He phase. 

After this brief introduction, the review is structured 
as follows: in §2 we will review the main theoretical dis- 
cussions and some experimental studies that appeared 
before the torsional oscillator (TO) "revolution" 
in §3 we will discuss the main experimental results in 
the last four years; section 4 is devoted to a discussion 
of the theoretical contributions (we will discuss mostly 
quantum simulation results) stimulated by the new ex- 
perimental findings, whereas §5 contains a discussion of 
what remains to be investigated and clariied on the basis 
of the present status of the research in this field. 

2. The Supersolid Phase before the TO "Revo- 
lution" 

2.1 The theoretical speculation 

ODLRO was discussed for the first time by Penrose 
and Onsager^) when they successfully generalized the no- 
tion of BEG for an interacting system of particles. The 
one-body density matrix pi{r,r') represents the proba- 
bility amplitude of destroying a particle in r* and creat- 
ing it in r. and ODLRO is defined by a nonzero value 
of pi in the limit |r — r*| oo. This implies a macro- 
scopic occupation of a single momentum state, i.e., BEG, 
because the Fourier transform of pi represents the mo- 
mentum distribution. In the same paper, by assuming a 
model for a solid in which particles are localized at their 
lattice sites (thus neglecting the possibility of exc;lianges 
and then explicitly breaking the Bose symmetry), they 
excluded the possibility of the simultaneous presence of 
crystalline order and ODLRO. This was not received as 
a definitive conclusion, and some years later the theoret- 
ical discussion concerned models for solids in which the 
assumption of the localization of the particles was weak- 
ened. The failure of the argument of Penrose and Onsager 
is only possible in a so-called quantum solid, i.e., a solid 
in which the mean-square displacement of the atoms is 
comparable to the interatomic distance. In a quantum 
solid the standard harmonic approximation completely 
breaks down, and it is no longer possible to associate an 
atom with a specific lattice position so that, even in the 
presence of a translational broken symmetry, the atoms 
are substantially dclocalized and microscopic processes 
such as exchange, the formation of vacancy-interstitial 
pairs, and more complex processes are allowed. 

Because of their low mass and weak interatomic inter- 
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action, ^He atoms are very difficult to solidify and, even 
at the lowest temperatures, they only solidify if a sub- 
stantial pressure is imposed (about 25 bar). Moreover, 
near the melting density the mean-square oscillation of a 
^He atom around the equilibrium position is a very large 
fraction of the lattice parameter, approximately 25%. 
Thus, solid ^He is the prototype of a quantum solid, and 
at the end of the 1960s it was readily recognized as the 
most likely candidate for observing a supersolid phase 
in nature. The mechanism for maintaining the atom de- 
localized and thus indistinguishable was proposed to be 
based on the presence of vacancies in the ground state of 
the solid phase: the so-called zero-point vacancies. This 
would mean that in the ground state the probability of 
occupation of a unit cell of the crystal is less than unity 
or, more precisely, two in the case of ^He, whic solidifies 
as a hep crystal at low pressure. Often a crystal in such 
a state is called an incommensurate crystal, in contrast 
with a commensurate crystal in which there is an inte- 
ger number of atoms in each unit cell of the crystal. If 
the atoms are bosons, then a vacancy, i.e., the absence 
of an atom, also obeys the same statistics, so that at low 
temperatures BEG can be expected with the presence 
of some form of superfluidity. All solids at a finite tem- 
perature contain a finite concentration of point defects 
such as vacancies, but this concentration in a classical 
solid vanishes exponentially fast as a function of 1/T, 
the inverse absolute temperature, as T approaches zero. 

Andreev and Lifshitz^^ suggested that vacancies, being 
mobile defects in solid ^He, behave as waves in the crys- 
tal, and have an allowed band of energies that correspond 
to different wavelengths (quasi-momenta) of these defec- 
tons; even if a localized vacancy costs a finite energy, the 
lower part of the band might be below the energy of the 
commensurate crystal. This happens if the width of the 
vacancy band is sufficently large, and this implies that 
the tunneling probability of an atom to the vacancy site 
is high. In this case, the ground state develops a finite 
concentration of zero-point vacancies, which can undergo 
BEG. This picture is based on the assumption of nonin- 
teracting vacancies, and Andreev and Lifshitz considered 
vacancy-vacancy interactions simply to limit their num- 
ber. As discussed below, the assumption of noninteract- 
ing vacancies appears to be incorrect, because even a few 
vacancies display significant mutual interaction. 

Ghcster's argumcnt^^ in favour of a supersolid phase 
of ^He was based on the demonstration by Reatto^^ that 
a certain class of wave functions, the Jastrow wave func- 
tions, always possess ODLRO, i.e., a finite condensate 
fraction. A Jastrow wave function for a system of A/" par- 
ticles in volume V is given by the following product of 
pair functions, 

*j(i?)=ne-^"('^-)/Qy' , (1) 

i<j 

where rij ~ |n — is the distance between particles i 
and j, R = {ri, ..,fv\/'} is the many-body configuration. 
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and Qj is the normahzation constant of ^fj, i.e., 




The function u(r), usuahy cahed the pseudopotential, in 
the description of "^He systems typically takes the McMil- 
lan form, u{r) = (b/r)"'', where b and m are positive vari- 
ational parameters. Such a simple form can be used to 
give a semiquantitative description of liquid "^He and, if 
the value of the "core" parameter b is sufficently large, 
describes a system in the solid phase. Such a class 
of wave functions has a configurational probability, i.e.. 
the square of the modulus of the wave function, which 
is equal to the configurational distribution function of a 
classical system at a fictitious temperature T* interact- 
ing via a two-body potential v*{r) proportional to the 
pseudopotential of ^j, i.e., v*{r) = kBT*u{r). As for 
any classical solid at a finite temperature, at a density 
where the system is solid such a configurational proba- 
bility must describe a solid with a finite concentration of 
vacancies;^) thus, for this model, zero-point vacancies are 
present and BEC is already proved.^' Some years ago, the 
concentration of vacancies in the state described by a Jas- 
trow wave function was estimated"'^^^ to be about 6x 10^^; 
however, one should keep in mind that gives a poor 
description of ^He in the solid phase, for instance, the 
mean-square deviation of an atom from the equilibrium 
position is about 1/3 of the correct value. A substantial 
improvement can be obtained by multiplying the Jastrow 
wave function by the product of one-body terms, usually 
Gaussians, which localize the atoms around equilibrium 
lattice positions, explicitly breaking the transla- 

tional symmetry; this is the Jastrow-Nosanow (JN) wave 
function: 

N 

^jN{R) = ^j{R)x^e-^\^^^-^^' /qy^ , (3) 

i 

where Qjn is the respective normalization. 

Following the works by Andreev and Lifshitz and by 
Chester, the presence of a supersolid phase was immedi- 
ately related to the possibility of measuring nonclassical 
rotational inertia by Leggett;-^'') a rotating soHd displays 
an inertia lower than that associated with the rigid ro- 
tation of all its atoms because a fraction of them do not 
follow the rotation of the container. Thus, the superfiuid 
fraction is given by Ps/p = (1 — I /Id), where Id is the 
classical moment of inertia, and the actual moment of 
inertia, 



is expressed as the limit of vanishing angular velocity w 
of the susceptibility of the free energy to the rotation. 
A rigorous upper bound for the superfiuid fraction was 
derived; it showed that at zero temperature the super- 
fluid fraction is reduced to a value less than 1 when the 
ground-state wave function lacks translational invariance 
and becomes zero in the absence of "exchanges" of par- 
ticles between lattice sites. Here a simple comment is 
worth mentioning because some successive contributions 



to the literature contained a misleading interpretation: 
being Leggett's discussion centered on an upper bound, 
the presence of "exchanges" does not provide a sufii- 
cient condition for the observation of NCRI. Leggett's 
letter also contained the suggestion of measuring NCRI 
by an experiment similar to the TO experiments which 
has been realized in 2004 and thereafter: the dawning of 
the revolution. 

2. 2 Lack of success in early experiments 

The interest in the possibility of a supersolid state in 
^He first motivated a number of theoreticians and then 
some physicists specialized in low-temperature experi- 
ments to study and search for quantum effects in solid 
^He. The research results obtained up to the beginning of 
the 1990s are reported in the review of Meisel.i°) During 
this period experiments did not reveal any signature of 
such a state, with the exception of some anomalies in ul- 
trasound measurements."'^*^ This led to a progressive loss 
of interest in this topic, and the supersolid state faded 
away as a theoretical dream. On the theoretical side, the 
approaches utilized in the 1970s were based either on 
very simplified models or on phenomenological models; 
thus, it was problematic to judge the relevance of their 
conclusions for solid ^He. The majority of these studies 
used quantum lattice models to investigate the proper- 
ties of the ground state and the excited states of quantum 
solids. Other approaches were based on the develop- 
ment of the procedure suggested by Leggett^^^ to obtain 
improved numerical estimates of an upper bound to the 
superfiuid fraction in the ground state of solid *He; these 
calculations used a phenomenological Gaussian model for 
the one-particle local density.^^"^''^ Some generalizations 
of the two-fiuid model were also attempted.^ ''■^^^ In the 
absence of experimental evidence for the supersolid state 
there was no strong motivation to apply the powerful 
simulation tools that were under development in those 
years to the search for the supersolid phase. 

The earlier experimental search for a supersolid phase 
was focused on the study of the plastic flow of objects 
moving in solid "^He.^^'^°) Later, the attempts to dis- 
cover such a state were based on mass flow^^'^^' or on 
TO experiments, ■^■'^ but no evidence of anomalies in the 
behaviour of solid ^He was found. The earlier TO experi- 
ment was carried out in the correct range of temperatures 
(T > 25 mK) and pressures (2-5 — 48 bar) but gave no 
indication of NCRI. However, under similar conditions, 
Kim and Chan were later able to observe NCRI; per- 
haps the sensitivity of the apparatus was inadequate or 
the annealing procedure reduced the signal to under the 
detectable level: it still was not time for a revolution. 

The only voice out of the chorus has been that of Good- 
kind and coworkers,^*' who reported that acoustic 
waves in solid *He are scattered by a nonphonon family of 
excitations with a small activation energy of the order of 
1 Some years later. Ho et al. observed a sharp peak 
in the acoustic attenuation in the presence of ^He impu- 
rities with a concentration of a few tens per 10^, and the 
anomaly decreased with increasing "^He concentration;-^^) 
these data were interpreted in a rather unusual way, in 
terms of a continuous second-order phase transition from 
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a Bose condensed state above a critical temperature of 
about 160—180 mK to a normal state below it. In another 
more recent experiment the interaction between trav- 
eling waves generated by heat pulses and sound waves in 
pure and impure hep solid ^He was investigated; the re- 
sults of this experiment were also interpreted as indicat- 
ing the presence of long-lived nonphonon excitations with 
a gap in their spectrum in the high-temperature phase, 
as previously suggested. It is worth mentioning these re- 
sults, even if they are still not understood, because the 
dependence on ^He concentration and the temperature 
range in which these phenomena appear seem to suggest 
a connection, if not a common origin, with the more re- 
cent quantum effects found in solid ^He. 

During those years, some attention was devoted to 
the measurement of the concentration of not only ther- 
mal vacancies, but also hypothetical zero-point vacan- 
cies. A summary of the measurements in solid helium is 
discussed in ref. 36. It was found difficult to reconcile 
the thermal vacancy content measured directly by X-ray 
diffraction with values based on the measurements of dif- 
ferent quantities (pressure, nuclear magnetic resonance 
of ^He impurities, ion mobility, plastic deformation, ul- 
trasonic attenuation). An output of these experimen- 
tal studies has been the vacancy activation energy; this 
quantity, measured by different methods near the melting 
density of solid ^He, was estimated to be in the range of 
7— 15 K, with the exception of the ultrasonic attenuation 
measurements of Goodkind^*^ (who assumed that the 
nonphonon family of excitations observed in his sound 
experiments is directly connected with vacancies), which 
gave an energy of about 1 K. In this context, zero-point 
vacancies have never been experimentally observed. In 
principle one can determine the presence of vacancies in a 
rather straightforward way. From the wave vectors of the 
Bragg peaks, one obtains the lattice parameters and the 
volume of the unit cell, which can be compared with the 
value obtained from the macroscopic density of the sys- 
tem. In reality, this is not so straightforward^'^^ because 
one never has a perfect single crystal completely filling 
a given volume. In addition, the direct X-ray method 
does not have the necessary sensitivity to measure a very 
small concentration. ■'^^ It is estimated that the present 
experiments cannot exclude ground-state vacancies at a 
concentration lower than about 0.4%.^^^ 

During this period, simulation techniques were used to 
study liquid and solid "^He from an accurate microscopic 
point of view, but without trying to answer the super- 
solid question. First attempts to obtain a quantitative 
description of solid ''He used variational techniques and 
the already introduced JN wave function.^^^ Among the 
variational techniques, it is worth mentioning a new class 
of wave functions, the shadow wave functions (SWFs), 
which were introduced in 1988 as a variational tool able 
to describe liquid and solid phases with the same func- 
tional form.^°) Using an SWF, correlations between par- 
ticles are introduced implicitly beyond the pair level of 
by coupling the particles to a set of subsidiary vari- 
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ables, called shadows, which are integrated over: 

^swF{R) = 'ij{R} [ dSG{R,S)^y{S)/Qy^p, (5) 
Jv 

where S = {si, .., sV} is the many-body configuration of 
the shadow variables, G{R,S) = J^f e-<^l^'-^'l' is the 
product of Gaussian one-body factors, '^j and '^j are 
Jastrow functions (in general with different pseudopo- 
tentials), which take into account the direct two-body 
correlations between particles {"^j) and between shadow 
variables (*j), and Qswf is the normalization 

QswF= [ dR [ dS [ dS'^'^j{R) x 
Jv Jv Jv 

xGiR,S)GiR,S')^'j{SmiS') . (6) 

An SWF is manifestly Bose symmetric, and the solid 
phase is described via a spontaneous broken translational 
symmetry process; when the density is above a certain 
value, the many-body correlations present in the SWF 
become so strong that the solid phase is stable and its 
energy is below that of the metastable liquid phase. It is 
important to note that this class of wave functions still 
provides the most accurate variational description of liq- 
uid and solid ''He, accurately predicting the freezing and 
melting densities;^^) the energy per particle obtained us- 
ing a fully optimized SWF is only 2 — 3% above that 
computed using "exact" QMC methods and with an al- 
most constant gap over a wide density range, giving rise 
to a very accurate equation of state. Chester's argument 
was also extended to this class of wave functions; SWFs 
thus describe a solid with a finite concentration of vacan- 
cies (which has recently been estimated^^^) and BEC.'*''^ 

"Exact" simulation techniques provided a more accu- 
rate approach to the quantitative microscopic study of 
the condensed phases of '^He. Results for the equation 
of state have been obtained at zero temperature by the 
Green's function Monte Carlo method^^^ (GFMC) and 
the diffusion Monte Carlo method'*i'45.46) (pMC). These 
zero-temperature methods are known to rely on the use 
of accurate variational trial states to implement the so- 
called "importance sampling", i.e., an improvement of 
the sampling algorithm, which restricts the exploration 
of the configurational space only to the "relevant" regions 
where the trial state does not vanish, thus substantially 
reducing the fluctuations in the computed quantum av- 
erages and allowing for their reliable calculation. This is 
a delicate aspect of the results obtained by these meth- 
ods because it not always easy to show that they are not 
affected by the choice of the trial state. 

This aspect is overcome in another "exact" QMC 
method, the path integral Monte Carlo (PIMC) method 
developed by Ceperley and Pollock (for a review see ref. 
47), which is suitable for computing finite- temperature 
quantum averages of a many-body system, 

(6) = Tr(p6)/Tr(p) , (7) 

where p = e~^^ is the density matrix. The method al- 
lows for the calculation of such averages by expressing 
the N-body density matrix in the coordinate representa- 
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tion at temperature T = {kBl3)~^ as a path integral: 

M-l M-1 

p{R,R',P)= / ]ldRi{l[{Ri\e-^"\Ri+,)} , (8) 

J i=l i=0 

where Rq = R and Rm = and by using an accurate 
approximation for the high-tcmpcratTire density matri- 
ces {Ri\e~'^^\R^_^_i) such that the quantum averages, 
computed via the Monte Carlo method, are not affected 
by this approximation within the statistical uncertainty 
of the calculation. For a boson system the symmetry of 
p{R, R', /3) is recovered via the direct sampling of permu- 
tations between single-particle paths, which correspond 
to ring polymers; not only the superfluid fraction, via the 
winding-number technique, but also off-diagonal proper- 
ties such as the one-body density matrix, which gives ac- 
cess to the Bose-Einstein condensate fraction, can also be 
computed. Finite-temperature properties have been 
studied by the PIMC method, which has been used 
to compute the kinetic energy''^) in liquid and solid ''He 
and also the Debye- Waller factor*^) in solid ^He with 
good agreement with experimental values deduced from 
neutron-scattering experiments. 

2.3 Early microscopic studies on a potential supersolid 

'^He phase 

Because of its ability to describe the solid phase via 

a spontaneous broken translational symmetry process, 
the SWF variational technique is particularly suitable 
to study disorder phenomena in quantum solids. For in- 
stance, the liquid-solid interface was studied"''"^ and the 
first microscopic quantitative estimation of the vacancy 
activation energy AEy was another application^^^ of the 
theory. AEy was foimd to depend strongly on the den- 
sity, and near the melting density AEy turned out to be 
about 15 K; both of these results were in agreement with 
the majority of the experimental estimations; the vacan- 
cies were also found to be mobile since they appeared at 
many different sites in the lattice during the simulation. 

Chester's conjecture on the relationship between va- 
cancies and BEC was quantitatively verified in 2001 us- 
ing an SWF in the presence of vacancies that were found 
able to induce BEC in the system. "''^^ The condensate 
fraction per vacancy, i.e., the condensate fraction divided 
by the concentration of vacancies = no/Xy, turned 
out to be about 0.21 ^He atoms per vacancy. In that 
work, different scenarios for the presence of a BEC phase 
in solid *He were discussed; four different phase diagrams 
were proposed depending on the presence of zero-point 
vacancies and on the values of quantities such as the va- 
cancy activation energy and the vacancy effective mass. 

Two years later, the same variational technique was 
used to compute the vacancy excitation spectrum and 
longitudinal phonon frequencies in solid ''He.^^'^'*^ The 
calculation was based on the development of the excited- 
state SWF variational technique used previously to study 
the phonon-maxon-roton spectrum in liquid ^He with 
very accurate results. The agreement of this microscopic 
and ab initio calculation of the vacancy spectrum with a 
tight-binding model was found to be quite good, and this 
suggested that a previously introduced hopping model 
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for a vacancy^^' can provide a reasonable approxima- 
tion when the correct bandwidth is known. From the 
computed vacancy excitation spectrum it was possible 
to extract the vacancy effective mass, which, in the hep 
lattice, turned out to have a small anisotropy and to be 
about 0.31 m4 in the basal plane and 0.38 m4 perpendic- 
ular to it, where TO4 is the mass of one *He atom. Never- 
theless, this value was a large amount of mass for nothing 
(a hole)! Also the bandwidth was found to be anisotropic 
in the range 7 — 11 K in hep solid "^He, depending on the 
direction of the wave vector in the crystal. The large 
value of the bandwidth A implies that atoms jump very 
frequently in the vacant site and that the residence time 
of a vacancy, estimated by r = h/A, is 0.6 x 10^^^ s, 
about 4 times larger than the period of a high-frequency 
phonon of the crystal. 

In 2003 the vacancy activation energy was com- 
puted by a projector quantum Monte Carlo method, 
the shadow path integral groimd state'''^"''®^ (SPIGS) 
method, developed from the path integral ground state 
(PIGS) method®®^ which is able to compute exact (within 
statistical errors) ground state (T = K) expectation 
values. The aim of the PIGS method is to obtain a state, 
^r, indistinguishable from the true ground state within 
the statistical uncertainty of the Monte Carlo calcula- 
tion by applying the imaginary time evolution operator 
e~^^ to a variational state, which in the SPIGS case is 
an SWF: 

^r = e-^"^swF . (9) 

By discretizing the path in imaginary time and exploiting 
the factorization property e^('^i+'^2)ff = p-riH ^-t2H ^ 
can be expressed in terms of convolution integrals that 
involve the imaginary time propagator {R\e~^'^^\R'), 
which, if St is sufficently small, can be approximated 
as is done with the PIMC method. This procedure maps 
the quantum system into a classical system of special in- 
teracting open polymers. If the starting variational state 
is Bose symmetric, as in the case of an SWF, the PIGS 
method automatically preserves this symmetry; thus, the 
calculations obtained by the SPIGS method in the solid 
phase of ^He provided the first "exact" ground-state cal- 
culations with the Bose symmetry preserved. Another 
appealing feature peculiar to the PIGS method is that 
in the variational ansatz acts only as a starting point, 
while the full path in imaginary time is governed by 
^-StH ^ which depends only on H. Therefore, the role 
of the starting variational wave function in the PIGS 
method is drastically reduced in comparison with other 
"exact" r = K methods such as GFMC and DMC. 
Vacancy activation energies computed by the "exact" 
SPIGS method were found in good agreement with the 
values obtained previously by the SWF technique.'"'^' 

This was the theoretical and experimental research 
scenario for solid ^He before the TO "revolution" . 

3. Recent Experimental Findings: the Revolu- 
tion Starts and Continues 

In the discussion of more recent experimental results 
on the investigation of a supersolid phase of ''He, we will 
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discuss experiments not in chronological order but by 
genre in order to provide a more complete view of the 
findings. We start with the TO experiments. 

3.1 Torsional oscillator experiments 

The scenario was radically changed by the paper by 
Kim and Chan in 2004,^^^ in which they reported the 
measurement of the oscillation period of a TO containing 
solid helium confined in Vycor glass, which is a porous 
material permeated by interconnected channels with a 
typical diameter of 60 A. "'He inside Vycor glass at low 
pressure and temperature is a superfiuid liquid. Under 
compression, ^He solidifies at a pressure larger than that 
required for a bulk sample. Kim and Chan performed the 
TO measurement in Vycor glass, and they found that be- 
low a temperature of about 0.2 K the period tq of the os- 
cillator dropped below the value at a higher temperature. 
Since tq = \/I /G /27r, where / is the moment of inertia of 
the TO and G is the torsion spring constant, if G remains 
constant this result implies that the moment of inertia 
decreased below that corresponding to the rigid rotation 
of the system. The fact that there was a critical-velocity 
eff'cct (at large amplitudes of the oscillator there was no 
missing inertia) and that the effect was suppressed when 
the ^Ho boson atoms were replaced by ^He fermions, led 
to the interpretation of these results as a manifestation 
of NCRI, sought for many years, associated with the su- 
persolid state in which about 1% of the ^He atoms did 
not respond to the imposed oscillation of the container, 
i.e., the superfiuid fraction PsI p was approximately 1%. 
The critical velocity turned out to be particularly small, 
of the order of few /um per second and of order of the ve- 
locity associated with one quantum of vorticity, and also 
this seemed to fit in with the supersolid scenario. These 
results was announced in 2003 during a talk by Chan 
at the Quantum Fluids and Solids International Sym- 
posimn (QFS2003) in Albuquerque (New Mexico, USA), 
news that magnetized the interest of many theoretical 
and experimental physicists present at the conference. 
The overall impression was that the confinement of solid 
**He was the key ingredient for inducing the necessary 
degree of disorder in the lattice, probably in the form of 
vacancies, to lead the system to a supersolid phase;. Only 
a few months later this scenario was cast into serious 
doubt by a second paper by Kim and Chan;^^^ the same 
effect was measured in bulk solid '*He. In this study the 
authors also noted that the effect was essentially sup- 
pressed by the presence of a plug in an annular con- 
tainer; the revolution had definitively started. It seemed 
that the two measures (Vycor and bulk) should possess 
a common mechanism at the origin of NCRI. In 2005 a 
similar effect was reported by the same group for solid 
**He confined in porous gold.^^-* It is important to note 
that the surface-to-volume ratio available for solid ^He 
changes by many orders of magnitude in these experi- 
ments; the presence of a similar phenomenon in the same 
range of temperatures under such diff'erent confinement 
conditions is an observation that still needs a reasonable 
explanation. This missing inertia is a robust effect not 
limited to the lowest pressure region where the solid is 
stable but occurs up to the highest pressures in present 
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experiments of approximately 140 bar.^^^ A puzzling re- 
sult is that the NCRI fraction in bulk "'He as function of 
pressure has a maximum at a pressure above melting, at 
about 50 bar.^-"-) In the light of the wide variation of the 
NCRI fraction measured later in different high-quality 
crystal samples (see below), these results should not be 
so surprising and might be due to an increased amount 
of disorder as pressure is increased. 

The experimental results of Kim and Chan have been 
independently reproduced by other groups around the 
world. ^^"^^^ TO experiments have also been replicated 
under a number of different conditions in order to extract 
a greater amount of information.^^' ^^"^^^ At this point it 
is important to summarize the findings from the set of 
experimental results obtained by the TO technique. 

All the above-mentioned experiments report an NCRI 
fraction, i.e., the superfiuid density ps in the supersolid 
scenario, which vanishes continuously as the tempera- 
ture is increased, similarly to a smeared phase transition 
rather than a sharp one, as expected in a homogeneous 
system. In addition, there is a dissipation peak (i.e., a 
minimum in the oscillation amplitude of the TO), which 
appears with NCRI at a temperature where ps has the 
maximum rate of change. This is similar to what was 
measured in TO experiments on *He adsorbed on planar 
surfaces. ''^^ The quality of the crystalline state of solid 
^He strongly depends on the growth method. Growth by 
the blocked-capillary method is known to produce poly- 
crystalline samples of varying sizes and defect contents, 
depending on the conditions of the growth and the an- 
nealing procedure. Growth at a constant pressure or at 
constant temperature from the superfiuid phase is known 
to produce single crystals or almost so. In all but one 
TO experiments, the solid was grown by the blocked- 
capillary method, thus the solid sample was polycrys- 
talline. 

Annealing a sample of solid *He is a standard pro- 
cedure to improve the quality of the crystal, and the 
effects of annealing have been reported by Rittner and 
Reppy.''^'^^'^*^^ They have eben been able to reduce the 
NCRI to below the level detectable in their experiment. 
The infiuence of annealing on NCRI has been also ob- 
served in other experiments;^^' "^^^ sometimes the overall 
effect has been to increase the NCRI fraction, and some- 
times it has reduced it. However, in contrast with the 
findings of Rittner and Reppy, it was never reduced to 
under the level detectable in the respective experiments. 
This indicates that disorder plays an important role in 
the phenomena. In fact, a large range of NCRI fractions 
(from 0.04% to about 20%) have been observed, with 
the larger values being observed in samples that presum- 
ably contain more defects. However, the amount of dis- 
order present in the samples was not determined quan- 
titatively in these experiments. This disorder is likely to 
be due to grain boundaries, vacancies, dislocations, or 
stacking faults. In one experiment, solid samples grown 
at constant pressure or temperature were used, and also 
in these cases the NCRI signal was detected, although 
with a reduced amplitude. On the other hand, high- 
quality crystals, all exhibiting NCRI, also gave consid- 
erably different values for the NCRI fraction depend- 
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ing on the sample. As mentioned above, the NCR! ef- 
fect also does not set in sharply at a precise "transi- 
tion" temperature, but it continuously increases up to 
a constant value with decreasing temperature; in high- 
quality crystals^^) the transition is much sharper and the 
NCRI fraction was found to be in agreement with the 
expected two-thirds power law for a continuous normal- 
superfluid transition (see Fig.l on J. Phys. Soc. Jpn., 
Vol.77, No. 11, p. 111010). However, as can be seen in 
Fig.l (on J. Phys. Soc. Jpn., Vol.77, No.ll, p.lllOlO) 
also in the high-quality crystals there is a residual knee 
at the temperature where NCRI sets in. Recently, re- 
markable hysteretic behaviour has been found in TO ex- 
periments,^^) with the NCRI phenomenon depending on 
the history of the oscillation velocity and on thermal cy- 
cling from a temperature above the transition tempera- 
ture to a low temperature. Another study of the thermal 
histories of the resonant frequency of a TO filled with 
solid ^He resulting from changes in the oscillation speed 
above and below the NCRI onset temperature has also 
been published. It showed that shifts in the resonant 
frequency are irreversibly induced by the variation of the 
oscillation speed below the onset temperature, revealing 
the presence of metastable states; however, the temper- 
ature dependence of the frequency shift for each sample 
is fully reproducible during warming and cooling scans 
at low oscillation speeds. Trapped vorticity might be at 
the origin of such metastability. 

Natural He contains about 0.3 part per million (ppm) 
^He impurities. NCRI is extremely sensitive to the con- 
centration of '^He. Additional '^He is found to broaden 
the onset of NCRI and to shift it to a higher "tran- 
sition" temperature but above a ^He concentration of 
order 0.1% NCRI disappears. More recently ultrapure 
samples of "^He have been studied with concentrations 
of ^He in the range 1 — 300 parts per billion (ppb). At 
such small concentrations the effect of '^Hc is also observ- 
able on the NCRI signal (see Fig.l on J. Phys. Soc. Jpn., 
Vol.77, No.ll, p. 111010). A reasonable explanation for 
this phenomenon is that the "'He atoms arc trapped by 
the defects of the crystal. For instance, there is evidence 
from ultrasound and shear modulus measurements''^' 
that a ■'He atom is bound to the dislocation core by an 
energy in the range 0.7-3 K. 

Rittner and Reppy^^^ have also observed an interest- 
ing correlation between the minimum-maximum NCRI 
measured in different TO experiments: it increases with 
increasing surface-to-volume ratio, S/V, which is specific 
to each TO cell. This has been interpreted as an indica- 
tion of the role of disorder in NCRI because a cell with a 
greater S/V, being more confining, should stabilize more 
defects in the crystal. 

Using a double-resonance TO^^'^^-' Aoki et al. were 
able to show that the reduction of the NCRI fraction is 
characterized by a critical- velocity effect and not by the 
displacement amplitude or the acceleration, because in 
these latter cases the reduction does not coincide in the 
two modes; moreover, they found that the observed dissi- 
pations and frequency shifts of the two modes cannot be 
made compatible with a modeF^^ in which NCRI is ex- 
plained in terms of the transformation of a liquid compo- 
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nent into a glass component at lower temperatures. This 
model also does not explain the different behaviour ob- 
served in the TO experiments in the presence of a blocked 
annulus. 

The value of the NCRI fraction varies by more than 3 
orders of magnitude in different samples. Also the range 
of temperatures in which NCRI has been detected varies 
but to a much smaller extent, i.e., in the range 60 — 
200 mK. When an experiment is conducted on ultrapure 
high-quality crystal, this range appears to be invariably 
reduced to 60 — 80 mK. 

Since the early measurements were performed in highly 
polycrystalline solids, grain boundaries have been sus- 
pected to play an important role in NCRI. The fact 
that crystals grown at a constant temperature or pres- 
gyj-g69) exhibited NCRI means that grain boundaries 
are not essential for NCRI, and other defects must be 
important. This also excludes the possibility that the 
metastable glassy region of "^He, which may be present in 
a rapidly quenched sample, '^^^ can play the decisive role 
in these phenomena. Dislocations are reasonable candi- 
dates, since it is known that even in a single crystal the 
density of dislocations can vary by a large amount de- 
pending on experimental conditions. A joint experiment 
on TO and on the content of dislocations in a single 
crystal is highly desirable. This major role of defects in 
the NCRI of solid ^He has even led to the conviction 
that the whole issue of the supersolidity of ^He is exclu- 
sively due to the presence of extrinsic defects, which any 
solid sample contains. By growing solid *He in the mK 
region, one can obtain single crystals with a very small 
number of dislocations, '^''^ at least dislocations ending on 
the crystal surface. It will be very important to perform 
TO measurements on such crystals. In this respect it is 
important to observe wether high-quality solid ^He sam- 
ples pass the bio eked- annulus test with a reduction of 
the NCRI fraction compared with that measured in an 
open cell, as theoretically expected for a supersolid. 

3.2 Elastic a.nom.alies in solid '^He 

Recently, a new anomaly has been found in solid 
^He.''*) The elastic shear modulus has been measured 
and a stiffening of the solid has been observed at low 
temperatures at which NCRI occurs, with an increase 
of the shear modulus which is approximately 10%. The 
T dependence of the shear modulus has a remarkable 
similarity with that of NCRI; also the "^He concentra- 
tion plays a significant role on the shear modulus, and 
has an effect completely similar to what it is observed 
for NCRI in TO experiments (see Fig. 2 on J. Phys. Soc. 
Jpn., Vol.77, No.ll, p. 111010). It has been proposed that 
this anomaly is due to dislocations, ''^^ i.e., the same kind 
of disorder that has also been discussed in relation to 
NCRI,®^' ''^' with dislocations being immobile at low tem- 
peratures, due to the pinning on "'He impurities, and be- 
coming mobile at higher TJ^^ On the other hand, should 
this completely explain NCRI, it would seem counterin- 
tuitive; the ability to move without dissipation seems at 
odds with the stiffening of the solid. The discovery of 
this stiffening is an important new information on the 
system, but we are far from understanding its cause and 
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its relation, if any, with NCRI. 

3.3 Mass flow experiments 

Liquid ^He in the superfluid state has a number of 
pecuUar properties in addition to NCRI. Obviously, in 
solid ^He other signatures of superfluidity have also been 
looked for. One method has been to perform experiments 
in which the measurement of mass flux inside the solid is 
attempted. In the first experiment of this kind, Greywall 
did not find any indication of siiperflow upon applying 
a small pressure difference between two chambers joined 
by a capillary and filled with solid ^He.^^^ Some years 
later, a second experiment in which a pressure difference 
was induced via different liquid-solid interface levels in a 
U-tube32) gave negative results. 

The observation of dc flow in solid ^He would provide a 
strong confirmation of superfiuid effects; therefore, soon 
after the first positive TO experiments, searches for such 
a fiow restarted. In 2005 Day et al.^^^ searched for the 
presence of a pressure-induced flow of solid ''He conflned 
in Vycor glass pores; neither this experiment nor the fol- 
lowing one,*"^ which was devoted to studying the flow 
of solid ''He through an array of capillaries, were able 
to detect any signature of such a flow. The main char- 
acteristic of these two experiments is that the pressure 
was directly applied on the solid lattice, thus raising the 
question of wether the correct "component" was being 
pushed. 

In 2006 a different version of the U-tube experiment 
succeeded in producing observable frictionless mass flow 
through the solid;^^) this phenomenon was observed only 
with the contemporaneous presence of some grain bound- 
aries directly detected in the cryostat at the liquid-solid 
interface. The range of temperatures in which this ef- 
fect was measured also extends above 1 K suggesting the 
possibility that this effect, even if very interesting, could 
not be related to NCRI in TO experiments. The original 
interpretation of this result was in terms of the superflu- 
idity of grain boundaries being the origin of the effect; 
however, the authors have recently proposed a different 
mechanism: the flow of mass along liquid channels 
present along the contacts between grain boundaries and 
walls, thus regarding the effect as a superflow of the liq- 
uid. 

Very recently, the results of a new experiment have 
been announced:^''^ Ray and Hallock have observed un- 
usual mass transport in hep solid "^He with an apparatus 
that allows the injection of *He atoms from the super- 
fluid directly into the solid using two Vycor glass rods, 
which enter the cell filled with solid ''He, and exploiting 
the fact that *He inside Vycor glass crystallizes at a much 
higher pressure than in the bulk. The flow was observed 
both in crystals grown by the blocked-capillary method 
and in crystals grown at constant temperature. This flow 
turns out to be pressure- and temperature-dependent; no 
flow is observed above about 27 bar and above 800 mK. 
At low pressure and at T below about 400 mK the flow 
induced by the unbalanced pressure between the two Vy- 
cor glass rods is approximately linear in time, as expected 
for a superfluid moving at a critical velocity. This char- 
acteristic temperature of 400 mK is at least two times 



higher than the characteristic NCRI temperature in the 
case of low-quality crystals and four times that for high- 
quality crystals. In this experiment there seems to be no 
space for explanations different from a mass flow of ^He 
atoms inside the bulk solid. For this reason, should these 
results be confirmed, this experiment may definitively 
show the; prescincc; of superflow in solid ^Hc. However, 
given also the different region of characteristic temper- 
ature and pressure, any relation with NCRI seen in TO 
experiments is far from being proven and understood, 
and further flow measurement at lower temperatures is 
urgently needed. 

3.4. Thermodynamic properties 

If the establishment of NCRI represents a true phase 
transition, one should expect to find some signature of 
it in equilibrium thermodynamic properties such as, for 
example, the specific heat. Early solid ^He heat capacity 
measurements hardly reached a low temperature limit 
comparable to the NCRI temperature range; moreover, 
the sample cells used in those experiments were con- 
structed with materials (e.g., heavy wall metal) that gave 
a large contribution to the heat capacity at low temper- 
ature.^^' This was also the case in the experiment by 
Clark and Chan,*^' where the authors extended the spe- 
cific heat measurement down to 80 mK without finding 
any sharp features. More recently, a peak in the specific 
heat has been reported at the temperature of the NCRI 
phenomenon^^' tising a sample cell made of silicon, which 
has a heat capacity more than 10 times lower than that 
of solid ''He in the considered temperature range (see 
Fig.3 on J. Phys. Soc. Jpn., Vol.77, No.ll, p.lllOlO). 
This is a very difficult experiment, and it will be impor- 
tant to obtain an independent confirmation of this find- 
ing. Recently, more accurate measurements have been 
performed; they confirm the presence of this peak, and 
its amplitude depends on the '"'He concentration being 
smaller in ultrapure samples of solid ''He (1 ppb). The 
presence of different concentrations of ^He does not alter 
the position of the peak, which is aromid 70 — 80 niK. The 
specific heat peak also depends on the growth method, 
with solids grown by the blocked-capillary method hav- 
ing a stronger peak. 

On the other hand, a very accurate measurement of 
the melting curve^'^"^^' in the temperature range 10 — 320 
mK found no contribution to the solid entropy other than 
the usual contribution due to phonons. The slope of the 
melting curve is proportional to the entropy difference 
between liquid and solid at coexistence, which at very low 
temperatures should be dominated by phonons. There- 
fore, the melting curve gives information on the excess 
entropic contribution in one of the two phases, and it 
should give sign of a phase transition, being the entropy- 
variations connected to the specific heat as function of 
the temperature. In the first experiment^''-' an anomaly 
in the melting curve below 80 mK was observed, but in 
a later analysis the authors claimed that it was entirely 
due to an anomaly in the elastic modulus of Be-Cu, which 
their pressure gauge was made of.^*' They also found a 
contribution to the entropy proportional to T"^ at low 
temperatures.^^' Such a contribution is at the limit of 
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their accuracy and the nature of this term in tlic en- 
tropy of solid ^He is still unclear. Anderson ct al., on 
the basis of a phenomenological theory, have proposed 
that such an entropy term may be due to strongly corre- 
lated vacancies. Recently, Maris and Balibar^^) have 
observed that this term can also be explained by anhar- 
monic phonon effects. 

The absence of an anomaly in the melting line is in se- 
rious disagreement with the presence of a peak in the spe- 
cific heat since both quantities probe the entropy of the 
system, and no explanation for the disagreement has yet 
been given. Excluding experimental artefacts, one possi- 
bility is that the anomaly in the entropy is absent in the 
difi'erence between the entropy of the liquid and that of 
the solid, since the melting line depends on this differ- 
ence. If this is the case we do not know its origin. 

3.5 More unexplained anomalies 

Hints of an elusive comprehension of solid ^He have 
also been obtained from inelastic neutron-scattering ex- 
periments. In 2002, neutron-scattering measurements on 
bcc solid ^He showed the presence of new unexplained 
nonphonon "optic- like" excitation modes, ^^'^^-^ and in 
2006, macroscopic structural fluctuations in bcc solid ^He 
were observed via neutron diffraction studies.^'*'' While 
a study of the Debye- Waller factor in hep solid ^He 
has shown the absence of any low-temperature anoma- 
lies down to 140 mK,^^^ recently, anomalous excitations 
have also been measured by inelastic neutron-scattering 
experiments in hep solid ^He.^^^ The authors have iden- 
tified a branch of clearly observable excitations with a 
quadratic dispersion curve centered on a reciprocal lat- 
tice point: a roton-like mode in solid *IIe, which was iden- 
tified as delocalized vacancy modes. In this respect, a the- 
oretical study^'') of ground-state commensurate ^He has 
shown the presence of vacancy-interstitial pairs as a fiuc- 
tuation effect. In variational SWF theory such pairs are 
unbounded, whereas such pairs turn out to be bounded 
in the exact SPIGS method. Since the ground state rep- 
resents the "vacuum" of any excitation of the system,*'^) 
it has been suggested that such vacancy-interstitial pairs 
represent a signature of a new kind of excitation different 
from phonons. 

It is also interesting to remember a study on the free 
expansion of solid ^He into vacuum^^' ^^^^ from a hole of 
diameter in the range 1 — 5 /ini in the container. First, a 
sequence of periodic intensity bursts was observed with 
a very regular period, which increases with increasing 
pressure or decreasing temperature, the so-called geyser 
effect. Above the lambda point at freezing, 1.76 K, the 
periodic bursts were observed only on solidification. At 
temperatures below 1.76 K, the pressure at which the 
geyser effect could be observed was 4-5 bar above the 
melting curve. Closer to the melting curve the flow is 
continuous in time and pure Bernoulli liquid flow is ob- 
served for a pressure of about 2 bar above the melting 
curve, suggesting that the solid closest to the melting 
curve possesses different "flow" properties. Low con- 
centrations of ^He, from 1% down to 0.1%, are sufficient 
to completely remove the "anomalous" behaviour near 
the melting curve, and the geyser effect is always observ- 



able at the melting curve down to the lowest measured 
temperature, about 1.3 K.^°°) In this experiment the sys- 
tem is far from equilibrium, and excess vacancies were 
proposed by the authors as a possible mechanism for the 
geyser effect, but even if these vacancies could play a role 
in the geyser effect, the different flow properties close to 
the melting line remain unexplained. We believe that it 
will be important to extend these experiments to a lower 
temperature, in the region of NCRI. 

4. Microscopic Studies after the TO "Revolu- 
tion" 

An ensemble of ^He atoms represents a strongly in- 
teracting many-particle system, which cannot be studied 
analytically via a perturbative approach with the aim of 
obtaining a quantitative microscopic picture of its physi- 
cal properties. The variational method, at least at T = 
K, is a viable method but the boson statistic of the atoms 
also makes the many-body problem exactly solvable, at 
least for an equilibrium state, both at T = K and at 
flnite T, by robust simulation methods (quantum Monte 
Carlo, QMC), even if one has to be careful because the 
results refer to a finite number of particles, usually less 
than 1000, with periodic boundary conditions. There- 
fore, the results of QMC mathods are always affected by 
finite size effects that should be removed before making 
a comparison with experiments. Depending on the quan- 
tity, the size effect can be more or less important, and 
it can be estimated by performing computations for sys- 
tems of different size. If the system is in an ordered solid 
state, there is another size effect: given the periodicity 
of the crystalline state and the use of periodic bound- 
ary conditions, there is a commensuration effect between 
the two periodicities, and this is much more diflacult to 
disentangle. 

QMC methods therefore play a leading role in the mi- 
croscopic theoretical study of the condensed phases of a 
strongly interacting many-particle quantum system such 
as liquid and solid *He. QMC methods have been exten- 
sively used in the past to acquire valuable information, 
especially on the properties of the superfluid phase of 
^He. Since the TO "revolution", all this expertise has 
been focused on the solid phase. A number of studies 
have also appeared in the last few years on solid ^He 
based on simplified models and phenomenological ap- 
proaches. In the present review we will mainly discuss 
theoretical results obtained using QMC techniques. 

4^.1 Quantum simulations: commensurate solid 

PIMC methods have been applied to study the off- 
diagonal properties of solid ^He at finite temperatures, 
and the result of such computations is that a commen- 
surate solid has a one-body density matrix, which de- 
cays exponentially at large distances; thus, the commen- 
surate crystal does not possess ODLRO.^°^' ^"^^ Also ps 
turns out to be zero. The fact that these results, even if 
not computed at temperatures below 0.2 K, do not show 
trends in the range 0.2 —IK exclude the possibility of a 
phase transition at lower temperatures in the proximity 
of T = 0.2 K. The PIMC method has also been used also 
to calculate exchange frequencies in bulk hep *He;^°^' 
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Fig. 1. (Color online) One-body density matrix computed in com- 
mensurate solid ^He. Different colours represent different algo- 
rithms used: (red squares) points taken from ref. 11 obtained 
by the "worm algorithm" PIMC technique, (blue lines) SPIGS 
method for different t values, (green lines) PIGS method applied 
to a JN wave function for different t values. (Upper figure) 2D 
triangular lattice at p = 0.0765 A"'^. (Lower figure) 3D hep solid 
''He near melting. 

by fitting the computed frequencies to a lattice model it 
was concluded that the commensurate crystal is not su- 
perfluid. 

Even if an SWF describes a quantum system with a 
finite concentration of vacancies, by choosing the num- 
ber of atoms to fit the periodic boundary conditions, one 
can also use this technique to simulate a commensurate 
crystal, and ODLRO has been found with a very small 
(5 X 10^®) condensate fraction. This discrepancy be- 
tween the "exact" PIMC method and the variational re- 
sults has been recently explained by computing pi by the 
"exact" SPIGS method. The most extensive compu- 
tations have been performed in a two-dimensional (2D) 
triangular lattice. Here the projector e~'^^ in eq. (9) was 
applied to two completely different wave functions. The 
first was an SWF, which is translationally invariant and 
has BEC, and the second was a JN function, eq. (3), 
which is not translationally invariant, has no BEC, and 
is not Bose symmetric. Note that the application of the 
Bose symmetric projector e~'^^ projects ^jat onto the 
subspace of Bose symmetric wave functions. As shown in 
upper Fig. 1 both computations give the same one-body 
density matrix pi within statistical errors, and piir^r') 
turns out to have dominant exponential decay with in- 
creasing |r — 7^1 so that there is no BEC. In the three di- 
mensional case (3D) (see lower Fig. 1), pi obtained start- 
ing form SWF and from JN show an exponential decay, 
similar to what was found with PIMC. Beyond approxi- 
mately 20 A , pi obtained starting from SWF and from 
JN slightly differ, and this suggests that convergence is 
essentially reached, but larger r values are needed in or- 
der to obtain complete convergence in this range, a ma- 
jor computational task. pi{r — f') in Fig. 1, computed by 
the SPIGS and PIGS methods, is shown for c? = r — 
along the nearest-neighbour direction, whereas pi com- 



puted at T = 0.2 K by the PIMC "worm algorithm" is a 
radial average, which explains the reduced oscilla- 
tions superimposed on the exponential decay. 

We conclude that there is overall agreement on the 
fact that commensurate solid ''He does not possess BEC. 
In a very recent computation by the DMC algorithm^°^^ 
a vanishing superfluid signal was found, but an extrap- 
olated estimator for the one-body density matrix gave 
a finite condensate fraction for the commensurate solid. 
This seems to contradict the result that ODLRO implies 
superfluidity;^' perhaps the discrepancy is due to the use 
of an extrapolated estimator for pi. It was also shown 
also that a vacancy is able to induce superfluidity, but 
the effect of correlations between vacancies has not been 
studied. 

Exact ground-state one-body local densities computed 
by the SPIGS method has been used to compute upper 
bounds for the superfluid fraction using an optimized 
one-body phase function. ^°^' The upper bound for com- 
mensurate solid ^He at melting was significantly higher 
(about 20%) than that experimentally observed, suggest- 
ing an improvement to the theory by permitting a many- 
body phase function; a theory for a two-body phase func- 
tion has already been developed. ^"^^ 

An interesting theoretical problem is why an SWF, a 
wave function that gives extremely accurate results for 
a number of properties of solid ^He, gives qualitatively 
incorrect results for the long-range limit of pi . It has been 
suggested^^^-* that the present SWF misses some long- 
range correlations, specifically those associated with the 
zero-point motion of transverse phonons. 

4-.2 Quantum simulations: vacancies and other defects 
The calculation of the one-body density matrix, pi, 
in the presence of vacancies has been obtained by the 
SPIGS projector method.^o^) It was found that a vacancy 
in a small fee crystal (107 atoms, 108 lattice sites) with 
periodic boundary conditions is able to induce BEC, and 
that the case of two vacancies was marginally compatible 
with a condensate fraction proportional to the vacancy 
concentration; these results are very similar to those ob- 
tained previously by the SWF technique. However the 
system was too small to be conclusive, and multiple va- 
cancies in the hep crystal were later found by PIMC 
method to be strongly correlated, such to form a 
bound state and phase separation has been proposed. ^^"^ 
Also computations by the SPIGS method"^) at T = K 
show that vacancies are strongly interacting but they do 
not form a single compact cluster for up to at least 6 va- 
cancies. The presence of attractive correlations between 
vacancies is also partially supported by an estimate of 
the elastic interaction potential between two static va- 
cancies in hep solid ''He; on the basis of an elasticity 
theory, an attractive interaction has been found along 
the basal plane, but the interaction is repulsive in some 
different directions. ^'^''•' 

Recent calculations of the vacancy effective 
mass^"'^' '^^^^ are in reasonable agreement with the 
early direct variational SWF calculation. Also the 
vacancy excitation spectrum^^^^ seems comparable. 
Glass states, grain boundaries, and a screw disloca- 
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tion have been studied in ^Hc systems by an imple- 
mentation of the PIMC method in the grand canonical 
ensemble.^^^'^^^'^^^^ The glass states were obtained by 
equilibrating the simulation at high temperatures and 
then "quenching" the system down to T = 0.2 K. This 
procedure leaves the system in a metastable amorphous 
"solid" state which possesses ODLRO (up to 0.5%) and 
a large superfluid fraction (up to 0.6). Concluding that 
the stochastic dynamic of a Monte Carlo sampling algo- 
rithm can resemble the real dynamic of a physical system 
exceeds the limits of what is reasonable; however, these 
results provide a useful insight of what metastable dis- 
order may be able to induce in solid ^He. Grain bound- 
aries have been found to be generically superfluid, with 
the exception of highly symmetrical ones; the transition 
temperature to the superfluid state has been estimated 
to be of the order of 0.5 K, thus about 10 times greater 
than the "transition temperature" of NCRI in TO exper- 
iments on high-quality crystals. Also, the simulation of a 
screw dislocation has recently revealed the capability of 
such a defect to sustain a superfluid flow along its core; 
however, in order to sustain an NCRI signal of about 1%, 
the authors^ have estimated on the basis of the results 
of their simulation that the required dislocation density 
should bo about 10^^ cm~^, many orders of magnitude 
greater than the observed experimental range (see for 
example ref. 117). 

Very simplified models of ''He confined in porous glass 
have been studied by the PIMC method^^®^ and the vari- 
ational SWF technique. 12°) The PIMC results sug- 
gest the presence of a persistent liquid layer just above 
the highly localized first layer of '*He atoms on the glass 
wall; this disordered layer gives rise to a superfluid re- 
sponse. Variational calculations were performed in a nar- 
row pore with radius i? = 13 A, and gave evidence that 
by increasing the number of atoms in the pore, solidifi- 
cation takes place layer by layer starting from the pore 
wall. Furthermore, the disorder trapped in the system is 
able to induce BEC even under complete solidification. 

The overall scenario of solid *Hc originating from the 
microscopic quantum simulation approach seems com- 
patible with a system which, if commensurate, does not 
possess the intrinsic property of superflow, but it is suf- 
ficient to include nearly any kind of disorder to obtain a 
system that can support some type of superflow. We can 
say that the commcnsiiratc solid is a marginal "insula- 
tor" , in the sense of the absence of BEC and superflow. 
The study of solid ^He systems by quantum simulation 
methods has given a useful insight into the effect of dis- 
order, but simulation of realistic disorder at zero tem- 
perature or at temperatures of the order of 50 mK will 
require the study of extended systems and considerable 
computational resources. 

4.-3 Commensurate or ineommensurate ground state 

Vacancy activation energies have also been computed 
by a grand-canonical PIMC technique^^") and more re- 
cently by the standard PIMC method. '^-''^) These calcu- 
lations substantially agree with previous variational^^-* 
and exact ground-state calculations^^' (the difference 
of about 3 K at melting is probably due to the different 
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He- He potential employed); thus, all microscopic QMC 
methods agree on the fact that the introduction of an 
extra vacancy in the ground state costs a finite amount 
of energy greater than 10 K. Moreover, these results do 
not show any particular dependence on the system size 
for up to about 500-1000 ^He atoms.^^' 1^°) This result 
of a finite energy cost for a vacancy has been taken as 
evidence that no vacancies can be present at low T, but 
the interpretation of these calculations is far from being 
unanimous."' 9"' 109'"°) 

The nature of the ground state, commensurate or in- 
commensurate, is a very important point in the discus- 
sion on a possible supersolid phase in solid "^He. We con- 
cluded that the microscopic calculations on solid ^He in 
literature do not allow one to infer wether the ground 
state is commensurate or incommensurate. i"^-' Obviously 
this has no more relevance for the Andreev-Lifshitz- 
Chester (ALC) scenario^'^^ of essentially independent 
zero-point vacancies because of the strong correlations 
found between these defects; however, it may have rel- 
evance for the existence of an incommensurate ground 
state with correlated defects; thus, in the following we 
again propose our considerations. 

Even if the ground state has zero-point vacancies, 
adding one additional vacancy is expected to increase 
the energy; this is similar to the case when in superfluid 
^He one adds one extra phonon to the ground state that 
contains the zero-point motion of phonons.^®-* Therefore, 
we expect that the excitation energy of the system con- 
tains a vacancy band with a gap, but this by itself is not 
sufficient to exclude an incommensurate ground state. In 
addition, in a QMC simulation of a crystalline solid, there 
is a commensuration effect as mentioned above. The rel- 
evance of such an effect was shown by the following ex- 
periment. ^^^^ Take a cubic volume and 108 ^He atoms at 
a density where the solid is stable; the exact ground state 
is a commensurate fee crystal. Now if one places 107 ^He 
atoms in the same volume, the exact ground state is in- 
commensurate; one observes that the number of lattice 
sites remains 108, one more than the number of particles. 

All QMC methods exploit the intrinsic statistical na- 
ture of quantum mechanics by computing quantum av- 
erages via the Monte Carlo sampling of a niuiiber of con- 
figurations of the system that have a nonvanishing prob- 
ability, as measured by the configurational probability 
function, which corresponds to the square of the mod- 
ulus of a wave function at zero temperature or to the 
diagonal part of the density matrix at a finite temper- 
ature. In order to diminish the commensuration effect 
one should consider a macroscopic system. One can ar- 
gue that all the microscopic approaches that describe the 
solid phase in a natural way via a spontaneously broken 
translational symmetry (Jastrow or SWF for the varia- 
tional method, SPIGS or PIMC for the "exact" methods) 
have a configurational probability function that gives a 
nonzero probability, not only for configurations that cor- 
respond to a commensurate crystal but also for those 
that correspond to an incommensurate crystal. This is 
due to the fact that a vanishing probability is only ob- 
tained for configurations with two or more particles at 
a distance lower than the hard-core radius of the inter- 
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action potential. Defects such as vacancies, interstitials, 
and also grain boundaries, dislocations and more com- 
plex defects in the solid phase do not give rise to such 
configurations with overlapping atoms, so that such de- 
fects have a nonvanishing probability of being present in 
the ground state. 

Thus, configurations corresponding to a commensurate 
crystal as well as configurations corresponding to an in- 
commensurate crystal will contribute to the normaliza- 
tion of a Jastrow wave function Qj (cq. (2)), of a shadow 
wave function Qswf (eq. (6)) hut also to the normaliza- 
tion of the SPIGS wave function ^f^, and to the bosonic 
partition function of ^He in the PIMC 

Zb = ^J2 [ dRp{R,rR,P) (10) 

(V in eq. (10) represents the sum over all possible TV! 
permutations imposed by the Bose statistic). The point 
is: which is the region of the configurational space which 
gives the dominant contribution to the normalization for 
the macroscopic system we are considering? For Qj and 
Qswf, Chester's argument proves that the overwhelm- 
ing contribution originates from regions of the configu- 
rational space where the crystal is incommensurate, and 
this is also true for Zb in eq. (10) for any finite T. Re- 
cently, we have computed the equilibrium concentration 
of vacancies in the ground state of solid ^Hc described 
in terms of an SWF.^^^ Close to the melting density we 
obtained X„ = 1.4 x 10~^, a value much larger than that 
found for the Jastrow case. Our value is below the present 
experimental upper bound, but it is quantitatively signif- 
icant. For instance, should the vacancies organize them- 
selves in dislocations, one would find a dislocation den- 
sity of no — 10^° cm^^ when X„ ~ 10^'^. The value 
of X„ obtained from Qswf is rather large but not \m- 
reasonable; however, one should be cautious. First, using 
the SWF, vacancies are not strongly interacting^^' as in 
exact computations. Second, wc do not know what hap- 
pens to given by the SPIGS ^'r in the limit r oo. 
Similarly we do not know what happens to Zb in eq. 
(10), in the T ^ limit. Wc conclude that at present 
the question of whether the ground state of solid '*He 
is commensurate or incommensurate is still unanswered 
and that answering this question is a priority. 

4- 4- Phenomenological theories 

In two recent papers Anderson^^^' suggested that 
the phenomenon observed up to now is not supersolidity, 
which might occur at a lower temperature, but a dissipa- 

tive regime dominated by quantized vortices. An analogy 
is made with high-temperature superconductivity above 
the superconductor phase transition in a regime where a 
Cooper-pair gap is still present even if the phase fluctu- 
ations destroy standard superfluidity. In this regime the 
system is better described in terms of an incompressible 
vortex fluid, which manifests itself, for instance, with an 
anomalous diamagnetic susceptibility, that is the quan- 
tity for a charged system which is the analogous of NCRI 
in the neutral system. Some aspects of NCRI in solid ^He, 
such as the critical velocity and some observed metasta- 



bility phenomena, are consistent with such a vortex fluid 
view of this state. In this picture, the system should be 
considered^^^) as a density wave rather than a solid, with 
the helium atoms flowing as a fluid through this density 
wave. 

On the basis of a phenomenological Landau theory, 
with the assumption that the hypothetical supcrsolid 
phase transition is continuous, Dorsey et al.^^'^' stud- 
ied the transition in *He by assuming that superfluidity 
is coupled with the elasticity of the crystalline lattice. 
They found that this assumption does not afi^ect the uni- 
versal properties of the superfluid transition and that a 
lambda-anomaly in the specific heat should appear. They 
also predicted anomalies in the elastic moduli near the 
transition and local variations in the superfluid transition 
temperature due to inhomogeneous lattice strains. 

Toner has recently used Landau theory to study the 
effect of quenched edge dislocations on a possible super- 
solid transition. ^^^' He found that quenched dislocations 
can promote the transition from a normal solid to a 
supersolid even if the dislocation-free crystal is normal 
down to T = K. 

On the basis of a Ginzburg-Landau theory, Ye^^^' has 
studied the possible phase transitions in ^He and the 
conditions for the existence of a supersolid phase; more 
recently he has investigated the elementary low-energy 
excitations inside a supersolid. 

The quantum mean fleld approximation for bosons at 
T = K leads to the Gross-Pitaevskii (CP) equation, 
a nonlinear Schrodinger equation. This equation is very 
successful in describing low-density cold atoms but has 
limited application to the strongly interacting superfluid 
"^He. Pomeau and Rica^^'') observed in 1994 that the GP 
equation with a nonlocal interaction kernel produces a 
crystalline state when the density is suflicently large. The 
dynamical properties of the model^^*' ^'^^^ display NCRI 
and a solidlike elastic response to a small stress. The 
nonlocality of the interatomic interaction has a key role 
in the model, but no account is given of the short-range 
correlations arising from the repulsive short-range inter- 
atomic interaction; thus, the results are interesting but 
it is unclear how relevant they are for solid ^He. 

4-5 Discussion: present understanding 

After the large number of experimental results pub- 
lished in the last four years, one can have confidence in 
the quantum nature of the anomalies measured in solid 
^He. However, what is the origin of the anomalies in solid 
^He? Are they intrinsic or extrinsic? Also, what is the na- 
ture of the ground state of solid ^He? Is it commensurate 
or incommensurate? The combination of the possible an- 
swers to these questions gives four possibilities; the fifth 
is the vortex liquid suggested by Anderson, but if the ex- 
periment by Hallock is confirmed, this case can probably 
be ruled out. Can we exclude any of the other possibil- 
ities? The case of an incommensurate ground state that 
is supersolid only due to extrinsic disorder does not seem 
reasonable and we will not discuss it. 

We first consider, for solid **He, a commensurate 
ground state that possesses BEG or superflow as an in- 
trinsic property. As we have previously discussed, exact 
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microscopic simulations have convincingly shown that 
given the He-He interaction potential, a commensurate 
solid does not possess ODLRO; thus, we can exclude the 
intrinsic-commensurate case. Moreover, in refs. 130 and 
11 the authors propose a series of considerations, specific 
to solids in which the translational symmetry is sponta- 
neously broken, that support the conclusion that a com- 
mensurate supersolid has zero probability of being found 
in nature because it requires an accidental symmetry be- 
tween vacancies and intcrstitials, which arc the individ- 
uated necessary ingredients for a quantum solid to show 
superfluidity. We find such considerations vague and 
not very convincing. However, if we accept them, the 
ODLRO found in a commensurate solid described using 
an SWF^'^^ is consistent with this argument; in fact, the 
SWF state, which is the exact ground state of an un- 
known Hamiltonian (not that of a system of ''He atoms), 
has a finite concentration of vacancies, and its simulation 
in a commensurate state is equivalent to forcing it to pos- 
sess the previously discussed accidental symmetry, and, 
in fact, unbound vacancy- interstitial pairs are present in 
the SWF.''^) Also, the recently introduced symmetric JN 
wave function^°^^ cannot be considered a counterexam- 
ple; this state contains one-body terms that explicitly 
break the translational symmetry, and thus represents 
the ground state of a Hamiltonian that is not transla- 
tionally invariant. If the above-mentioned considerations 
are true, a system whose ground state corresponds to a 
commensurate solid cannot be supersolid; no arguments 
or counterexamples against this conclusion are known. 

We now consider the intrinsic-incommensurate case. 
This case not only includes the ALC scenario^' of 
zero-point vacancies, but also that in which the ground 
state contains a correlated state of vacancies^^^ or va- 
cancies and intcrstitials. One can also consider other 
defects^^^' or a correlated mixture of different defects. 
There is a consensus on the fact that tlu^ ALC sce- 
nario^' of essentially independent vacancies cannot 
alone be the explanation of the NCRI and of other 
anomalies recently observed in solid *He. Even if we ac- 
cept that the simple estimation of the vacancy activa- 
tion energy, giving a positive energy cost in systems with 
less than a few thousand atoms, does not exclude the 
presence of ground-state vacancies, ^"^^ the observation 
of strong correlations between vacancies with some indi- 
cation of a possible bound state^^""^^"*^ leaves only the 
possibility of a scenario in which, if the ground state con- 
tains vacancies, such defects are highly correlated. ^^^^ 
One can argue that this is not probable but there is no 
strong evidence to exclude it. The case of correlated va- 
cancies can be ruled out if phase separation occurs; our 
personal view is that present computations are not con- 
clusive regarding the issue of phase separation, and stud- 
ies of larger systems with increasing number of vacancies 
are needed. It is possible, for example, that vacancies are 
converted into other defects or that the continuous mu- 
tation into different defects occurs recursively. However, 
we should also consider the possibility that a hypotheti- 
cal supersolid phase might be associated not with vacan- 
cies but with some other correlated defects that, in any 
case, would make solid ^He different from the textbook 



picture in which each atom occupies a crystal cell. We 
know that the solid is periodic; this is demonstrated by 
the observation of Bragg peaks in the scattering of neu- 
trons or X-rays, but how the atoms populate the crystal 
cells may offer some surprise. The fact that we still do not 
have sufficient elements to exclude the incommensurate- 
correlated-defect ground state also indicates that we can- 
not use the microscopic results of quantum simulations, 
in which no ODLRO has been found for the commensu- 
rate case, to deduce that the ground state of solid ^He 
is commensurate and that the conditions of these simu- 
lations were not an accidental realization of the system. 

We now consider the remaining possibility, the 
extrinsic-commensurate case, i.e., the possibility that all 
the "supersolid" phenomena in solid *He (with a com- 
mensurate ground state) are due to the presence of some 
disorder induced by external conditions such as the con- 
dition of crystal growth. If this were the case, it would 
be of interest per se but, of course, not as interesting 
as if supersolidity were an intrinsic property. In such a 
case if a more and more perfect crystal is grown, NCRI 
would disappear and no effec;t will remain, leaving the 
ground state of solid ^He qualitatively similar to that of 
a normal solid. Microscopic quantum simulations have 
shown that almost any kind of disorder is able to induce 
superfluidity in the system; however, on the basis of the 
same results, no defect has been found to be capable of 
clarifying the observed anomalies. 

As mentioned above, the evidence from exact simula- 
tions is that commensurate solid *He has no ODLRO and 
no superfluidity. An interesting question is whether this 
result depends on the form of the interatomic interac- 
tion. It is possible that a commensurate solid might have 
ODLRO and superfluidity if the interatomic interaction 
differs substantially from that appropriate for "^He. 

5. Conclusions 

The recent surge of interest in solid ''He has provided 
much new information on the properties of this fascinat- 
ing solid. Anomalies now are not limited to TO evidence 
for NCRI but are also present in the shear modulus and 
speciflc heat. Solid *He has turned out to be much more 
diiEcult and interesting than expected, and we do not 
yet have a clear overall picture of the properties of this 
system. All attempts to explain the behaviour of solid 
*He in more conventional terms have failed so far, and 
the presence of a nonconventional state appears to be the 
most likely explanation. Should the report of superflow 
in solid ^He be confirmed, this conclusion will be defini- 
tive. However, the nature of this nonconventional state is 
still unclear. Experiments as well as microscopic theories 
indicates that the TO results, even if they are due to su- 
persolidity, cannot be interpreted in terms of the original 
ALC scenario of weakly interacting vacancies. Vacancies 
are strongly interacting and experiments clearly indicate 
that other defects, such as dislocations, grain boundaries, 
and ^He impurities, must have an important role. Actu- 
ally, the measured quantities depend so sensitively on the 
growth method, annealing, and ^He content that some 
researchers believe that all the anomalies are only due to 
extrinsic defects. The question of whether such a noncon- 
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ventional state is indeed only due to extrinsic defects or 
whether it is an intrinsic property of this quantum sohd 
is a fundamental question that requires an answer. An 
argument in favour of the intrinsic nature of the state 
is as follows. Although, as already mentioned, the value 
of the NCRI fraction varies by many orders of magni- 
tude depending on the experimental details, the tran- 
sition temperature is affected only marginally by such 
details. If supersolidity were entirely due to extrinsic de- 
fects, it would be rather difficult to understand why the 
transition temperature is hardly affected when the num- 
ber of defects varies by many orders of magnitude. In 
fact, the presence of quantum effects also in high-quality 
^He solid samples points in the direction of a noncon- 
ventional state of strongly correlated-defects; a model of 
this kind, based on a network of dislocations was devel- 
oped in ref. 124, and the transition temperature turned 
out to depend strongly on the defect concentration. Ob- 
viously, this argument does not exclude the possibility 
that some other extrinsic-defect model might agree with 
experiments; we simply observe that a more natural ex- 
planation comes from the other possibility as follows. If 
the underlying ground state has BEC, the system will 
have its own transition temperature but the intrinsic 
value of the NCRI fraction may be very small. The pres- 
ence of extrinsic defects might have an amplifying effect 
on the NCRI fraction without necessarily markedly af- 
fecting the transition temperature, i.e., extrinsic defects 
might change the NCRI fraction without significantly af- 
fecting other properties of the system. The present the- 
oretical understanding is that a commensurate crystal is 
normal even at zero temperature and that this commen- 
surate normal state is marginally stable, i.e., almost any 
deviation from it leads to a superfluid response. There- 
fore, if the anomalies of solid ^He below a temperature 
in the range 70 — 200 mK are due to intrinsic supersolid- 
ity modified by the presence of extrinsic disorder, ^He in 
the ground state ground state cannot be a commensurate 
crystal. Some researchers argue; against this possibility 
because, according to present computations there is an 
energy cost upon introducing any of the defects consid- 
ered so far. Others researcher argue that this is not con- 
vincing because one is neglecting the commensuration ef- 
fects present in quantum simulations. In our opinion the 
question of whether the ground state is commensurate or 
incommensurate is still open as a fundamental question, 
and hopefully experiments with improved sensitivity will 
give an answer to whether intrinsic disorder exists even 
at the lowest temperatures. From theory, we also expect 
extended investigations in a more realistic way of defects 
like multiple vacancies or dislocations as well as the study 
of the microscopic dynamics of the system. 

From the experimental point of view, it is important 
to perform joint measurements on the same sample us- 
ing different techniques in order to characterize solid "^He 
more completely. The expertise of growing crystalline 
^He of the highest quality should be used to investigate 
supersolidity phenomena. These are just a few examples 
of studies that should be undertaken. The issue of super- 
solidity will probably remain at the centre of attention 
in future. We have learnt something but much more has 



to come. 

Note added in proof: During the 2008 summer new 
experimental results have been announced; in the follow- 
ing, we report on some of them. Rittner and Reppy have 
performed an enhanced version of the bio eked- annulus 
TO experiment which confirm that NCRI effects are due 
to inertia and not to stiffness variations of the solid sam- 
ple (A. S. C. Rittner, workshop "Supersolid 2008", 18- 
22 August 2008, ICTP, Trieste, Italy). As reported by 
Beamish (J. Beamish, workshop "Supersolid 2008", 18- 
22 August 2008, ICTP, Trieste, Italy), the shear modulus 
anomaly is not at the origin of NCRI and is not an ef- 
fect of statistics because the anomaly is present also in 
hep ^He and it is essentially the same as in hep ^He. On 
the other hand, no anomaly is present in bee ^He. TO 
experiments (group of M. H. W. Chan) on solid ^He sam- 
ples, grown from the same helium gases used by Beamish, 
show no NCRI signal from either hep and bee solid ^He. 
Also, more accurate specific heat measurements in solid 
^He with different ^He concentrations indicate that the 
magnitude of the specific heat peak does not depend on 
the ^He concentration but it does depend on the crystal 
quality (M. H. W. Chan, ULT 2008, 14-17 August 2008, 
Royal HoUoway University of London, London, UK). 
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